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Abstract 



We construct the complete SU(2) parity-violating (PV) it, N, A interaction 
Lagrangian with one derivative, and calculate the chiral corrections to the 
PV Yukawa NNit coupling constant through 0(1/A^) in the leading 
order of heavy baryon expansion. We discuss the relationship between the 
renormalized h n , the measured value of h w , and the corresponding quantity 
calculated microscopically from the Standard Model four-quark PV inter- 
action. 



PACS Indices: 21.30.+y, 13.40.Ks, 13.88.+e, 11.30.Er 



I. INTRODUCTION 



The parity-violating (PV) nucleon-nucleon interaction has been a subject of interest 
in nuclear and particle physics for some time. To date, PV observables generated by 
this interaction remain the only experimental windows on the AS = 0, nonleptonic weak 
interaction. Since the 1970's, the PV NN interaction has been studied in a variety of 
processes, including p—p and p— nucleus scattering, 7-decays of light nuclei, the scattering 
of epithermal neutrons from heavy nuclei, and atomic PV (for a review, see Refs. |]||. The 
on-going interest in the subject has spawned new PV experiments in few-body systems, 
including high-energy p—p scattering at COSY, n+p — > at LANSCE [||], — > n+p 
at JLab ||, and the rotation of polarized neutrons in helium at NIST. 

The theoretical analysis of these PV observables is complicated by the short range 
of the low-energy weak interaction. The Compton wavelength of the weak gauge bosons 
(~ 0.002 fm) implies that direct and Z exchange between nucleons is highly sup- 
pressed by the short-range repulsive core of the strong NN interaction. In the conventional 
framework, longer range PV effects arise from the exchange of light mesons between nucle- 
ons. One requires the exchange of the tt, p and ui in order to saturate the seven spin-isospin 
channels associated with the quantum numbers of the underlying four-quark strangeness- 
conserving PV interaction, Tt^(AS = 0) (henceforth, the AS = will be understood). 
These exchanges are parameterized by PV meson-nucleon couplings, h M , whose values 
may be extracted from experiment. At present, there appear to be discrepancies between 
the values extracted from different experiments. In particular, the values of the isovector 
7riViV coupling, h n , and the isoscalar pNN coupling, h° - as extracted from p — p scatter- 
ing and the 7-decay of 18 F, do not appear to agree with the corresponding values implied 
by the anapole moment of 133 Cs as measured in atomic PV 0. 

The origin of this discrepancy is not understood. One possibility is that the use of p 
and 10 exchange to describe the short-range part of the PV NN interaction is inadequate. 
An alternate approach, using effective field theory (EFT), involves an expansion of the 
short-range PV NN interaction in a series of four-nucleon contact interactions whose 
coefficients are a priori unknown but in principle could be determined from experiment. 
The use of p and u exchange amounts to adoption of a model - rather than the use 
of experiment - to determine the coefficients of the higher-derivative operators in this 
expansion. Whether or not the application of EFT to nuclear PV can yield a more self- 
consistent set of PV low-energy constants than the meson-exchange approach remains to 
be seen. A comprehensive analysis of nuclear PV observables using EFT has yet to be 
performed. 

The least ambiguous element - shared by both approaches - involves the long-range 
7r-exchange interaction. At leading order in the derivative expansion, the PV ttNN in- 
teraction is a purely isovector, Yukawa interaction. The strength of this interaction is 
characterized by the same constant - h n - in both the EFT and meson-exchange ap- 
proaches. At the level of the Standard Model (SM), is particularly sensitive to the 
neutral current component of H.™ ' . In this respect, the result of 18 F PV 7-decay mea- 
surement is puzzling: 
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K = (0.73 ± 2.3K 



(1) 



where g n = 3.8 x 10~ 8 gives the scale of the h M in the absence of neutral currents |J. This 
result is especially significant, since the relevant two-body nuclear parity- mixing matrix 
element can be obtained by isospin symmetry from the /3-decay of 18 Ne 0. The result 
in Eq. ([[]) is, thus, relatively insensitive to the nuclear model. Theoretical calculations of 
h v starting from Ti™ have been performed using SU(6) tt symmetry and the quark model 
|7],§] , the Skyrme model || , and QCD sum rules |I| . As a benchmark for comparison with 
experiment, we refer the SU(6)^/quark model analysis of Refs. J7],[S| — hereafter referred 
to as DDH,FCDH []. These authors quote a "best value" and "reasonable range" for the 

/^(best) = 7g w (2) 
K (range) : (0 -> 30)^ . (3) 

where here the "best value" is more aptly described as an educated guess, while the 
"reasonable range" indicates a set of numbers such that theory would be very hard-pressed 
to explain were the experimental value not found to be within this band. Neverthess, the 
difference between the "best value" of Eq. (fj) and the 18 F result would appear to call 
for an explanation, and in the following note we comment on a possible source of the 
discrepancy. 

In general, the problem of relating the fundamental weak quark-quark interaction to 
the low-energy constants which parameterize hadronic matrix elements of that interaction 
is non-trivial. In the framework of EFT, one may define these constants at tree-level in the 
hadronic effective theory. The quantities extracted from experiment in the conventional 
analysis, however, are not the tree-level parameters, but rather renormalized couplings. 
Denoting the latter as h EFF , one has 

h EFF = z N Jz~ n h BARE + Ah 7T , (4) 

where h BARE is the coefficient of the leading-order, PV Yukawa interaction in the effec- 
tive theory, \[Z$ and \[Z K denotes chiral loop renormalizations of the nucleon and pion 
wavefunctions, respectively, and Ah n denotes contributions from chiral loops and higher- 
dimension operators to the Yukawa interactions (only the finite parts of these couplings 
are implied; loop divergences are cancelled by the corresponding pole terms in h BARE and 
the Zn^). At leading order in 1/A X , one has Zn >7T = 1, Ah n = 0, and h EFF = h BARE . The 
renormalized coupling appears as the coefficient in the one-pion-exchange (OPE) PV NN 
potential 



*Note that although the DDH analysis used the symmetry group SU(6) W in order to connect 
weak vector meson and pion couplings the predictions relating pion couplings alone to hyperon 
decay data rely only SU(3). 
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where g NNlT is the strong 7rNN coupling and /^(r) = exp(— m 7r r)/47rr. Neglecting the 
effects of three-body PV forces and 27r-exchange interactions, it is h FFF to which the 
result in Eq. ([[]) corresponds. 

The relationship between h FFF and the coupling obtained by computing (NTr\l-t™\N) 
in a microscopic model is not immediately transparent. In what follows, we make several 
observations about this relationship. We first show that Z^y/Z^ and Ah^ are substan- 
tial, so that h FFF differs significantly from h^ ARE . To that end, we compute all of the 
chiral corrections to the PV Yukawa interaction through 0(1/A|), where A x = AttF^. We 
work to leading order in l/m N in heavy baryon chiral perturbation theory (HBChPT). Of 
particular significance is the dependence of Ah n on other low-energy constants parame- 
terizing PV 2tc production and the PV iV — * 7T7tA transition. We subsequently reexamine 
the SU(6) w /quark model calculation of Refs. [0,|J and argue that most - if not all - of the 
chiral loop effects which renormalize h n are not included in the microscopic calculation. 
Thus, the relationship between h FFF and microscopic calculations remains ambiguous at 
best. This ambiguity is unlikely to be resolved until an unquenched lattice QCD cal- 
culation of h n using light quarks becomes tenable. In the meantime, one should not 
necessarily view a discrepancy between the experimental value of h FFF and microscopic 
model calculations as disturbing. 

Our discussion of these observations is organized as follows. In Section 2 we sum- 
marize our conventions and notation, including the PV chiral Lagrangians relevant to h n 
renormalization. Section 3 gives a discussion of the loop calculations. In Section 4 we 
comment on the scale of the loop corrections and provide simple estimates of some of the 
new PV low-energy constants appearing in the analysis. Section 5 gives our discussion 
of the relationship between h FFF and the calculation of Refs. |7]f|. Section 6 summarizes 
our conclusions. Some technical details are relegated to the Appendices. 

II. NOTATIONS AND CONVENTIONS 



We follow standard HBChPT conventions 11,12 and introduce 



in 



e, e = exp(— ), n = -7r a r a 



(6) 



with F n = 92.4 MeV being the pion decay constant. The chiral vector and axial vector 
currents are given by 



(7) 
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For the A, we use the isospurion formalism [13|, treating the A field T^(x) as a vector 
spinor in both spin and isospin space with the constraint r*T^(x) = 0. The components 
of this field are 

n-m\ ^-(sn ■ w 

The field also satisfies the constraints for the ordinary Schwinger-Rarita spin-| field, 

7^T; = and p% = . (10) 

We eventually convert to the heavy baryon expansion, in which case the latter constraint 
becomes v = with the heavy baryon velocity. 

The relativistic parity-conserving (PC) Lagrangian for 7T, N, A interactions needed 
here is 

C c = ^-TrD^D^ + NiiV^ - m N )N + g A NA^ lB N 

+ |^^ 7 a 75 + |( 7 ^ + 4 7 ,)75 + f 7/ ^ 7 a 757,]77 
+97rNA[Tt(g IMU + ztfyflJwiN + h.c] , (11) 

where uo 1 ^ = tr[r*A Ae ]/2 while and T>^ are the gauge and chiral covariant derivatives, 
respectively. Explicit expressions for the fields and the transformation properties can be 
found in ||14||. Here, z$ is an off-shell parameter, which is not relevant in the present work 
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In order to obtain proper chiral counting for the nucleon, we employ the conventional 
heavy baryon expansion of C PC , and in order to cosistently include the A we follow the 
small scale expansion proposed in In this approach energy-momenta and the delta 
and nucleon mass difference S are both treated as small expansion parameters in chiral 
power counting. The leading order vertices in this framework can be obtained via P + TP + 
where Y is the original vertex in the relativistic Lagrangian and 

are projection operators for the large, small components of the Dirac wavefunction re- 
spectively. We collect some of the relevant terms below: 

C PC = N[iv ■ D + 2g A S ■ A]N - ifftiv ■ D lJ - 5 ij 5 + gi S ■ A ij ]T^ 

+g*NA[f?4N + Nu«T?] (13) 

where is the Pauli-Lubanski spin operator and 5 = dia — mjv- 
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The PV analog of Eq. flllD can be constructed using the chiral fields X£ R defined as 

IT 



X a L = ^T a £, X a R = ^, X a ± =X a L ±X a R . (14) 

We find it convenient to follow the convention in Ref. [13] and separate the PV Lagrangian 
into its various isospin components. 

The hadronic weak interaction has the form 

U w = ^=-Jx-J X 1 + h.c. , (15) 

where J\ denotes either a charged or neutral weak current built out of quarks. In the 
Standard Model, the strangeness conserving charged currents are pure isovector, whereas 
the neutral currents contain both isovector and isoscalar components. Consequently, Ti w 
contains AI = 0, 1, 2 pieces and these channels must all be accounted for in any realistic 
hadronic effective theory. 

We quote here the relativistic Lagrangians, but employ the heavy baryon projections, 
as described above, in computing loops. It is straightforward to obtain the corresponding 
heavy baryon Lagrangians from those listed below, so we do not list the PV heavy baryon 
terms below. For the irN sector we have 

£aLo = h Q v NA^N (16) 

£a/=i = ^NYNTr(A^Xl) - ^M^NTr^X*) (17) 
'^-.F^NX 3 N 



2^2 



AI2 = h 2 v l ab N[X a R A,X b R + X a L A,X b L }YN 
h 2 



-^N[X a R A,X R - X a L A^X b L ]n 5 N , 

where T ab is a matrix coupling the X a,b to I — 2, 13 — 0. The above Lagrangian was first 
given by Kaplan and Savage (KS) fT5| . However, the coefficients used in our work are 



slightly different from those of Ref. [T5| since our definition of A^ differs by an overall 
phase. 

The term proportional to contains no derivatives. At leading-order in 1/F n , it 
yields the PV NNn Yukawa coupling traditionally used in meson-exchange models for 
the PV NN interaction 0,0. Unlike the PV Yukawa interaction, the vector and axial 



vector terms in Eqs. (|16| - |18|) contain derivative interactions. The terms containing h\ 



and h\ start off with NNtttt interactions, while all the other terms start off as NNn. 
Such derivative interactions have not been included in conventional analyses of nuclear 
and hadronic PV experiments. Consequently, the experimental constraints on the low- 
energy constants h l v , h l A are unknown. 
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It is useful to list the first few terms obtained by expanding the Lagrangians in Eqs. 
( |T6[|T8| ) in 1/F n . For the present purposes, the following terms are needed: 

^Yukawa = -^(P™r + ~ np^l " 3^^" + \^)\ ( 19 ) 
Cf N = - h ° v ^ h2v iprnD,n+ + fvfpD^-] (20) 

TT 

h 1 - h 2 

+i A 2 A nYl 5 n{7r + D fl 7i- - tTIV + ) 

IT 

+i^P^n>K + Dy - i^p^ 5 n7: + Dy . (21) 

7T 7T 

For the PV ttNN Yukawa coupling we have also kept terms with three pions. 

The corresponding PV Lagrangians involving a iV — > A transition are somewhat more 
complicated. We relegate the complete expressions to Appendix A, and give here only 
the leading terms required for our calculation. As noted in Ref. [Q, the one-pion ttNA 
PV Lagrangian vanishes at leading order in the heavy baryon expansion. The two-pion 
terms are 

•7pA+ + 7r~7r° •7pA++7r°7r~ 

£\ NA = - A p2 pa;+d»k-*o - P A++D^\- 

n n 

■ipA + ir°TT -ipA + TT+TT- 



TT TT 

•jpA+7r _ 7T+ •TpA°7r+7r° 

A -pA+D»n-n + - pA° u D^ + n° 



p2 r A* p2 



■ipA°TT°TT + •> pA~7r + 7T + 

-pA°n^°7r + - pA-D'V-V 



P2 f M P2 

TT TT 

•/,nA ++ 7r _ 7r~ ■'LnA + -K~-K 

nAj+^Tr-Tr- - nA+^n-n* 

n A+D^ 7f- - ^— nA° D'W 

TT TT 

nA° D^tt+tt- - nA° D^7r-7r + 



7T 7T 

•/,nA~7r+7r° •/,nA _ 7r°7r+ 

J)A T) TT TT 

p 2 UL ^V ly 71 71 ^2 



^ nA;D^7T + 7r° nA~D^TT\ + + h.c. , (22) 



where the couplings h p ^ ++7T n are defined in terms of the various SU(2) PV low-energy 
constants in Appendix A. 
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The PV 7rAA Lagrangians, also listed in Appendix A, contain terms analogous to 
the Yukawa, V, and A terms in Eqs. (|I6|-|T8|). Since we compute corrections up to one- 
loop order only, and since the initial and final states are nucleons, the PV irirAA terms 
(A-type) are not relevant here. The leading, single- 7r Yukawa and V^-type interactions are 



Yukawa 



"A" 



7T 



7T 



hi 

V3 V 



-Z^(A°A-7T 



- A + A 

+ - A^A°7T 



2h 



A (A+ A°tt + - A u A 



o A+7r - 



AA 



fly 



(A 



ILy 



h 



V 



++ lfl A+D^ + - 
-(A+ 7m A />V 



A+ 7m A ++ /> 



7T 



(23) 



(24) 



where the coefficients are given in Appendix A. 

One may ask whether there exist additional PV effective interactions that could con- 
tribute at the order to which we work. In the pionic sector there exists one CP-conserving, 
PV Lagrangian: 

C = ^XW-^[) • (25) 

At leading order in 1/F n , C w contains five pions. Its lowest order contribution appears at 
two-loop order at best, so we do not consider it here. 

Similarly, one may consider possible contributions from two-derivative operators. 
There exists one CP-conserving, PV operator: 



A 



-Na^[D,A u 



(26) 



There exist three independent PC, two-derivative operators [ 16| . For example, one may 
choose the following three: 

^-Ni^D^N , (27) 



— NA»A U N 
A v 



(2f 



^Na^[A„A u ]N. 



(29) 



As we discuss in Appendix B, none of the two-derivative operators in Eqs. (|26|) - (|29]) 
contribute to the renormalization of at the order to which we work in the present 
analysis. 
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III. THE LOOP CORRECTIONS 



The leading order loop corrections to the Yukawa interaction of Eq. ( |I9|) are generated 
by the diagrams of Figs. 1-2. As we discuss in Appendix B, the contributions from many of 
the diagrams which nominally renormalize h n vanish at the order at which we truncate. In 
particular, none of the vector (V-type) ttNN and 7rAA terms contribute to this order. In 
what follows, we discuss only the non- vanishing Yukawa and A-type contributions. Details 
regarding the vanishing of the other contributions appear in Appendix B. Following the 
conventional practice, we regulate the loop integrals using dimensional regularization. 
The pole terms proportional to 1/D — 4 are cancelled by appropriate counterterms. We 
identify only the terms nonanalytic in quark masses with the loops. All other analytic 
terms are indistinguishable from finite parts of the corresponding counterterms. 

The nonvanishing contribution from Fig. 1 (a) arises from the insertion of the 3ir part 
of the Yukawa interaction of Eq. (|TT| ) . The nonanalytic term is 

*M (a) = §^]n(^) 2 ^T + , (30) 

where A x = A7rF n and fi is the subtraction scale introduced in dimensional regularization. 
For simplicity, we show here only the contributions for n — > pn^. The terms for p — > nir + 
are equal in magnitude and opposite in sign since it is the hermitian conjugate of the 
n —>■ piT~ piece. This property holds to all orders of chiral expansion. 

The nonvanishing contribution from Fig. 1 (b) arises from strong vertex correction 
to the leading order ttNN Yukawa interaction: 

^) = I^M^-/Kt+. (31) 

The terms in Figs. 1(d)- (c2) are generated by the PV axial imNN couplings proportional 
to the h A . We have 

3 

lM {cl)+{c2) = 2V21X g A ^j-h\r+ . (32) 

The contribution from h\ to these two diagrams cancels out, leaving only the dependence 
on h\. We note that although this term is propotional to m?J 2 and, thus, nominally 
suppressed, the coefficient of h\ is fortuitously large (~ 1/4). The two pion vertex in 
FIG. 1 (dl)-(d2) comes from the chiral connection V^. 

^M {dl)+id2) = -^M^yKT+. (33) 

The leading contribution involving A intermediate states arises from Fig. 2(a). The 
corresponding amplitude receives contributions from three different isospin combinations 
for the A intermediate states. Their sum reads 



8 



»M 2W = -H2%^[(2^ - m») In(iL)* - 4J V^^Ib ^ '"> + . (34) 
9 Af, m„- v m„- 

A. 

The corrections generated by the PV nirNA vertices are 



,2 



■H/r ^ 9-kNA u x2 3 2 2 2 2n3/2i ^ + V ^ "^ilA H 
?M 2(61)+2(62) = W ~ ^ ~ m -) 111 — ' 

' 1 X 



(35) 



where h\ is defined as 



(36) 

Summing all the non-vanishing loop contributions yields the following expression for 
Ah n : 

20 ^A [(252 _ <} ln( J^ )2 _ 45 /^T^ ln 6 + V 52 -^! 



9 LV w/ v m/ v w m 

X, 



.^9nNA ux 2 3 2 . /i 2 r>^2 2\3/2i ^ + V ^ "^irA /o^N 

+ 3^Af [(5 -2 mJ51n W ~ 2(5 ln ^ ^ (37) 

The final nonvanishing corrections arise from N and 7r wavefunction renormalization. 
These corrections, which have been computed previously fll7l , generate deviations from 
unity of and \J~Z~l appearing in the expression for h FFF in Eq. ([|). In the case of Z N , 
the nonvanishing contributions arise from Figs. I(el)-(e2) and 2(cl)-(c2): 

7 i _ 9 n 2 m l ]n( V ]2 An 2 ^5 2 -ml ii 2 



8\ 5 2 — ml 5 + \ 5 2 — m\ 
-4-*-— -In ± -] . (3* 

A2 m w 



The pion's wavefunction renormalization arises from Fig. 2 (k) [IH 



z,-i = -Ut£\ m(JL) . (39) 



3 V A x / \m T[ 



Numercially, the loop contributions to \[Z~l are small compared to those entering Z N . 
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Note the one loop renormalization of h n from the PV Yukawa nNN and 7rA A vertices 
is already at the order 1/A£. An additional loop will introduce a factor of 1/A£. Loops 
containing the axial vector NNirn and NAtttt vertices and one strong NNtc or NAtt 
vertex are of 0(1/ A^F^). To obtain contributions of 0(1/ A^), one would require the 
insertion of operators carrying explicit factors of 1/A X into one loop graphs. We find no 
such contributions. 



IV. THE SCALE OF LOOP CORRECTIONS 

We may estimate the numerical importance of the loop corrections to h BAHE by taking 
5 = 0.3 GeV, g A = 1.267 @ and g nNA = 1.05 @ and by choosing ^ = A x = 1.16 GeV§ 
With these inputs, the value of Z^y/Z^ is completely determined. The vertex corrections, 
which appear as Ah n in Eq. (f|), depend on the PV couplings h n , h\, h&, and h\. We 
obtain 

h EFF = 0.5h n + 0.25h l A - 0.24/ia + 0.079^ . (40) 

Note that the effect of the wavefunction renormalization corrections is to reduce the 
dependence on h BARE by roughly 50%. In addition, the dependence of h EFF on h\ and 
h& is non-negligible. Their coefficients are only a factor of two smaller than that of 
h BARE . Although these contributions arise at 0(p 2 , p 3 ), they are fortuitously enhanced 
numerically. Thus, in a complete anaysis of the OPE PV interaction one should not ignore 
these constants. 

At present, one has no direct experimental constraints on the parameters h\, h&, and 
h^, as a comprehensive analysis of hadronic PV data including the full chiral structure 
of the PV hadronic interaction has yet to be performed. Consequently, one must rely on 
theoretical input for guidance regarding the scale of the unknown constants. Estimates 
of h\ are given by the authors of Ref. fllSl . These authors observe that the usual pole 



dominance approximation for P-wave non-leptonic hyperon decays typically underpredicts 
the experimental amplitudes by a factor of two. The difference may be resolved by the 
inclusion of local, parity-conserving operators having structures analogous to the A-type 
terms in Eq. fll7|). The requisite size of the AS = 1 contact terms may imply a scale 
for the analogous AI = 1 PV terms. If so, one might conclude that h\ should be on the 
order of 10^. On the other hand, a simple factorization estimate leads to h\ ~ 0.2^. 
While the sign of h\ is fixed in the factorization approximation, the sign of the larger 
value is undetermined. Thus, it is reasonable to conclude that h\ may be large enough 
to significantly impact h EFF , though considerably more analysis is needed to yield a firm 
conclusion. 



^Since the dependence on fi is logarithmic, one may choose other values, such as /i 
without affecting the numerical results significantly 



m 



pi 
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The 7rAA Yukawa coupling Ha has been estimated in Ref. || using methods similar 
to those of Ref. |7|]. The authors quote a "best value" of h& = —20g n , with a "reasonable 
range" of (—51 — > 0) x g n ^ Naively, subsitution of the best value into Eq. (^) would 
increase the value of h FFF , whereas the 18 F result would seem to require a reduction. As 
we argue below, however, the relationship between the couplings computed in Refs. |7],§] 
and the parameters appearing in Eq. ppp is somewhat ambiguous. Direct substitution 
of the theoretical value into h FFF may not be entirely appropriate. 

To date, no theoretical estimate of the A-type 7T7tNA coupling has been performed. A 
simple estimate of the scale is readily obtained using the factorization approximation. To 
that end, we work with tree-level form of Ti^ . Neglecting short-distance QCD corrections 
and terms containing strange quarks, one has 

H P J{AS = 0) = % cos 2 9 c u lx (l - l5 )ddj x (l - l5 )u (41) 

-2(1 - 2sin 2 e w )V^A® x + ^sin 2 9 W V® > A® A } , (42) 

(3) (3) 

where V£ and A x denote the third components of the octet of vector and axial vector 
currents, respectively, and 

n (0) = l -{u lx u + d lx d) . (43) 



Consider now the first term in the expression for h\ given in Eq. fl3~6|). In the factorization 
approximation, Ti,^ contributes only to the antisymmetric combination 

The neutral current contribution to this combination, which arises only from the term 
containing V x 3 \ is 

V2G F F2(1 - 2sin 2 6 w )C£(nA°) « 2^C 5 A (nA°) , (45) 

where C^(nA°) ~ Oil) is the axial vector n — > A form factor at the photon point. After 
Fierz re-ordering, the charged current component of Ti^ contributes roughly 

- (4^/3)C 5 A (nA ), (46) 

yielding a total factorzation contribution of about (2g n /3)C^(nA°). Thus, one would 
expect the scale of the axial vector tcttNA couplings to be on the order of a few x g w . 

In the particular case of the combination appearing in h^, however, the sum of factor- 
ization contributions cancels identically. As one sees from the expressions for the h 
given in Appendix A, isospin requires 



NA 
A 



7V7Y 



^This coupling is denoted /a An in Ref. 
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The factorization contributions independently satisfy this sum rule. The second combi- 
nation of constants appearing in Eq. (|56]), 

^W.jjfW ( (48) 

also vanishes in the factorization approximation, even though the individual couplings 
do not. The third pair of couplings received no factorization contributions. Thus, one 
has h\ = in this approximation. In principle, non-factorization contributions yield a 
non-zero value for h^. Although we have not evaluated these contributions, we do not 
expect the scale to be significantly larger than the factorization value for the individual 
foNAnn COU pii n g S , Consequently, we estimate a reasonable range for of (0 — > few) x g n . 

These theoretical estimates suggest considerable ambiguity in the prediction for h FFF . 
In principle, some of this ambiguity might be removed by performing the comprehen- 
sive analysis of hadronic PV suggested above, in which the various constants would be 
determined entirely by experiment. The viability of such a program remains to be seen. 



V. COMPARING WITH MICROSCOPIC CALCULATIONS 



The results in Eqs. (|37|-^) embody the full SU(2) chiral structure at 0(p 3 ) of 
(N^ll-C^lN) at leading order in the pion momentum. Any microscopic calculation of 
this matrix element which respects the symmetries of QCD should display the depen- 
dence on light quark masses appearing in h FFF . In principle, an unquenched lattice QCD 
calculation with light quarks would manifest this chiral structure. In practice, however, 
unquenched calculations remain difficult, and even quenched calculations require the use 
of heavy quarks. For a lattice determination of (Nir\T-l™\N), the expressions in Eqs. 
( p7H39| ) could be used to extrapolate to the light quark limit, much as the chiral structure 
of baryon mass and magnetic moment can be used for similar extrapolations 



In the absence of a first principles QCD calculation, one must rely on symmetries 
and/or models to obtain the PV NNtt coupling. A variety of such approaches have been 
undertaken, including the SU(6) w /quark model calculation of Refs. [7||§, the Skyrme 
model |§, and QCD sum rules [jl(J. To date, the DDH/FCDH analysis remains the most 
comprehensive and has become the benchmark for comparison between experiment and 
theory. Consequently, we focus on this work study" in the problem of matching 

microscopic calculations onto hadronic effective theory. 

The DDH/FCDH approach relies heavily on symmetry methods in order to relate 
the PV AS" = matrix elements to experimental AS = 1 nonleptonic hyperon decay 
amplitudes. All the charged current (CC) contributions to the AS" = 0, 1 B — > B'M 
amplitudes, where M is a pseudoscalar meson, can be related using SU(3) arguments. 
Likewise, the neutral current (NC) component of the effective weak Hamiltonian belonging 
to the same multiplets as the CC components {i.e. those arising from a product of purely 
left-handed currents) can also be related via SU(3). The remaining NC contributions to 
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the AS* = PV amplitudes are computed using factorization and the MIT bag model. The 
DDH approach also employs SU(6) W symmetry arguments in order to calculate parity- 
violating vector meson couplings. Although one requires only SU(3) to determine the 
pseudoscalar couplings, we refer below to the general SU(6) W formalism used in Refs. 

SI- 

The general SU(6)^ analysis employed by DDH/FCDH introduces five reduced matrix 
elements: a t , v , b t)V1 and c v . These constants correspond to SU(6)„, components of the weak 
Hamiltonian: 

[(55) 35 ® M 35 ] 3 5 ~ c v (49) 

[(55) 405 ® M 35 ] 280> 280 ~ h, K (50) 
[(55)405 ® ^f35] 2 80,280 ~ a t, a v (51) 

One may represent these different components of Ti 1 ^ diagramatically as in Fig. 3. The 
components shown in Fig. 3a,b correspond to b tjV and c v , respectively. In practice, these 
contributions are determined entirely from empirical hyperon decay data. The term in 
Fig. 3a corresponds to a t)V and is computed in Refs. using factorization. 

The PV NNtt Yukawa coupling can be expressed in terms of these SU(6) W reduced 
matrix elements plus an additional factorization/quark model term. Temporarily neglect- 
ing short-distance QCD corrections to H™, one has 

1 

(pn \H W \n) = — ^=tan0 c c„ (52) 

2 1 

esc 26 l c sin 2 9 w (2c v — b t ) + -sin 2 9 w y , 



9V2 



where 8 C and 9 W are the Cabibbo and Weinberg angles, respectively, and y denotes a 
Fierz/factorization contribution. The first term on the RHS of Eq. fl52j) gives the CC 
contribution, while the remaining terms arise from weak NC's. Including short-distance 
QCD renormalization of Ti.^ leads to a modification of Eq. d52|) : 



(pn-\H™\n) = {[l-2sin 2 ^] 7 (ir) + sin 2 ^}-4-^ 
L > sin r 



where 



+ sin 2 # c (5 1 + 5 2 ) , (53) 



9n = —7= tan 9 c c v (54) 



3\/2 

Bi = -^VE(K) f . / ) (b,./C> b,/12 c,./2) (•>•)) 
9V2 Vsm0 c cos0 c ' 

B 2 = ±F(K)y , (56) 
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and j(K), E(K) and F(K) are summed leading log (renormalization group) factors de- 
pendent on 



a s (/x) ri1 2 Arll M, 



Jr= l__^ i[1 l__ W/]ln _i . (57) 

The overall scale factor p appearing in Eq. (|5^) was introduced in Ref. 0] in order to 
account for various theoretical uncertainties entering the analysis. 

The appearance of c v , b t , and b v in g n and B 1 relies on tree-level SU(6)„, symmetry- 
long-distance chiral corrections of the types shown in Fig. 4 have not been explicitly 
included. Inclusion of such corrections would necessitate a reanalysis of the AS = 1 
amplitudes in much the same way that one treats the octet of baryon axial vector currents 
fTT|j or magnetic moments [21|. For example, letting A(A°_) denote the amplitude for 



A — > pn one has at tree-level 

A(A°_) = ^=(b v /6-b t /12-c v /2) . (58) 

Including the leading chiral corrections would yield the modification 

1 



A(A°_) = - / =^Z A Z p Z 7r (b v /Q~b t /12~c v /2) + AA(A i i) , (59) 

where AA(A°_) denotes vertex corrections and possible contributions from higher- 
dimension operators. Similar corrections would appear in the SU(6) W symmetry terms in 
Eqs. ([52], |53| ). Given the absence of these corrections from the DDH/FCDH analysis, the 
symmetry components (p7T~\Ti. P ^\n) do not formally embody the subleading chiral struc- 
ture of h FFF . The numerical impact of applying chiral corrections to the DDH/FCDH 
tree-level SU(6) W analysis is much less clear, since some of the chiral modifications can be 
absorbed into renormalized values of the chiral couplings, which are determined empiri- 
cally. Nevertheless, the potentially sizeable effect of the SU(2) chiral corrections on h FFF 
should give one pause. 

A related issue is the degree to which ambiguities introduced by kaon and r] loops in 
SU(3) HBChPT could plague an analysis of the AS" = 1 amplitudes. Here recent work by 
Donoghue and Holstein argues that finite nucleon size call for long-distance regularization 



of such heavy meson loops, which substantially reduces their effects ||22|| . Results are then 
similar to what arises from use of a cloudy bag approach to such matrix elements ||23|| . 
A comprehensive study of such issues - and their impact on the DDH/FCDH calculation 
of h n - goes beyond the scope of the present work. Nevertheless, the sizeable impact of 
the chiral corrections in h FFF and the use of tree-level symmetry arguments in Refs. |7],§] 
points to a possibly significant mismatch between h FFF and h° DH . 

The remaining terms in the DDH/FCDH analysis - involving the parameters rj and 
y - are determined by explicit MIT bag model calculations. One may ask whether the 
latter effectively includes any part of the subleading chiral structure of h FFF . In order to 
address this question, we make three observations: 
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1. Sea quarks and gluons generate c v . The parameter c v vanishes identically in any 
quark model in which baryons consist solely of three constituent quarks. The AS = 1 
hyperon decay data, however, clearly implies that c v 7^ 0. In order to obtain a nonzero 
value in a quark model, one requires the presence of sea quarks and gluons. It is shown in 



24J] for example, that c v ^ when gluons are added to the MIT bag model. Similarly, one 
would expect contributions from the qq pairs in the sea. Since relativistic quark models 
already contain qq pairs in the form of "Z-graphs" p5| , it is likely that disconnected qq 



insertions (see Fig. 5b) give the dominant sea quark contribution to c v . In a chirally 
corrected analysis of nonleptonic decays, the long-distance part of the disconnected qq 
insertions appear explicitly in the guise of pseudoscalar loops, while the short-distance 
contributions are subsumed into the value of c v and possible higher dimension operators. 
"Quenched" quark models without explicit pionic degress of freedom generally do not 
include the long-distance physics of disconnected insertions. 

2. The m^-dependence is different. In conventional HBChPT analyses of hadronic 
observables, one only retains the loop contributions non-analytic in the light quark mass. 
The constituent quark model (without explicit pions) has a difficult time producing these 
non-analytic contributions. The simplest, illustrative example is the nucleon isovector 
charge radius, (r 2 )T=i, which is singular in the chiral limit f26| . This chiral singularity, 
of the form lnm^ ~ \nm q , is produced by n loops. Relativistic quark models, such as 
the MIT bag model, yield a finite value for {r 2 )x=i as m q — > 0. One cannot produce the 
chiral singularity in a quark model without including disconnected qq insertions dressed 
as mesons. 

The corresponding argument in the case of h FFF is less direct, but still straightforward. 
In the limit of a degenerate N and A, the non-analytic terms in h FFF have quark mass- 
dependences of the form m g \nm q or m^ 2 . As we show in Appendix C, bag model matrix 
elements of the four quark operators appearing in Ti™ have a Taylor series expansion 
about m q = 0. Thus, the parameters 7] and y cannot contain the non-analytic structures 
generated by the diagrams in Figs. 1-2. 

3. Graphs are missing. This observation is simply a diagrammatic summary of the 
previous two observations. For simplicity, consider a subset of the quark-level diagrams 
associated with the appearance of h l A in h FFF . Typical contributions to the axial NNtctt 
PV vertex are shown in Fig. 5a. The corresponding loop contributions to h FFF appear 
in Fig. 5b, c. Those in Fig. 5b involve disconnected qq insertions, which do not occur in 
the constituent quark model. The contribution of Fig. 5c involves Z-graphs, which are 
produced in a relativistic quark model^. In principle, the 3q + qq intermediate state could 
contain an Ntt pair. As argued previously, however, the Z-graphs implicit in the MIT bag 
model calculation of h n do not produce the nonanalytic structure of the corresponding 
7r loop. Apparently, only an unquenched quark model, which generates the disconnected 
insertions of Fig. 5b, could produce the requisite nonanalytic terms. 



^e.g., as a correction to the bt )V terms of Fig. 3b 
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From this "case study" of the DDH/FCDH calculation of h n , we conclude that the 
SU(6) w /quark model approach used in Refs. |7||| does not incorporate the chiral structure 
of h EFF . Were the numerical impact of the chiral corrections negligible, this observation 
would not be bothersome. The actual impact of the chiral corrections, however, is signif- 
icant. Thus, it is perhaps not surprising that the 18 F result and the DDH/FCDH "best 
value" do not agree. 

VI. CONCLUSIONS 

With the confirmation of the electroweak sector of the Standard Model at the 1% 
level or better in a variety of leptonic and semi-leptonic processes, one has little reason 
to doubt its validity in the purely hadronic domain. Similarly, the predictions of QCD 
in the perturbative regime have been confirmed with a high degree of confidence. Thus, 
one may justifiably consider TC^, the effective Hamiltonian including its perturbative 
strong interaction correction, to be well understood. Moreover, the precision available 
with present and future hadronic PV experiments is unlikely to match the levels achieved 
in leptonic and semileptonic processes. Consequently, one has little hope of detecting 
small deviations in TL™ from its SM structure due to "new physics" . On the other hand, 
much about QCD in the non-perturbative regime remains mysterious: the mechanism 
of confinement, the dynamics of chiral symmetry breaking, the role or sea quarks in the 
low-energy structure of the nucleon, and so forth. Each of these issues bears on one's 
understanding of matrix elements of 'H P Y . In this sense, the low-energy, PV hadronic 
weak interaction constitutes a probe of the dynamics of low-energy QCD, in a manner 
analogous to the probe provided by the electromagnetic interaction. 

From a phenomeno logical standpoint, the matrix element one may hope to extract 
from hadronic PV observables with the least ambiguity is (NttIT-C^IN). In this study, we 
have argued that any theoretical interpretation of this matrix element must take into ac- 
count the consequences of chiral symmetry. Indeed the chiral corrections to the tree-level, 
PV ttNN Yukawa coupling are not small. At (9(1/A^), the effective coupling measured in 
experiments, h EFF , depends not only on the bare coupling, h BARE , but also on new (and 
experimentally undetermined) PV low-energy constants, h\, h^, and h/±, as well. Fur- 
thermore, the coefficients of h BARE , h\, and h& are comparable in magnitude. At present, 
one has only simple theoretical estimates of the magnitudes of the h\ and in addi- 
tion to the FCDH calculation of h&. These estimates suggest that the new PV couplings 
appearing in h EFF could be as large as h BARE . Since no experimental constraints have 
been obtained for the new couplings, there exists considerable latitude in the theoretical 
expectation for h EFF . 

For two decades now, the benchmark theoretical calculation of (iV7r|7i^|iV) has been 
the SU(6) w /quark model approach of Ref. [0], updated in Ref. ||. We have argued, 
however, that the DDH/FCDH calculation does not manifest the general strictures of 
chiral invariance obtained in the present analysis. At the quark level, this chiral structure 
reflects the importance of the "disconnected" qq components of the sea. While relativistic 
quark models contain qq sea quark effects in the guise of Z-graphs or lower-component 
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wavefunctions, the most common "quenched" versions do not include explicit disconnected 
pairsQ. Given the size of the chiral corrections associated in part with the disconnected 
insertions, it may then be not so surprising to find a possible discrepancy between the 
experimental value for h FFF and the DDH/FCDH "best value". 

Applying chiral corrections to the SU(3) analysis of AS = 1 hyperon decays may 
help to close the gap between h FFF and h° DH . Presumably, similar corrections should be 
applied in other approaches not containing explicit pionic degress of freedom. In the longer 
run, one may be able to use the chiral structure of h FFF to extrapolate an unquenched 
lattice calculation with heavy quarks into the physical regime. 
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**Some effects of disconnected qq pairs may, however, hide in the effective parameters of the 
quark model, such as the string tension \27\. 
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APPENDIX A: PV LAGRANGIANS 



Here we present the full expressions for some of the PV Lagrangians not included in 
the main body of the paper. The analogues of Eqs. (|T^-[TB|) are 



£li=o = fie^Nt^XlA.Xl + X a R A,X R ]T» 

+ gi N[A^ X a _} + T» + g 2 N[A^ X a _}^ + h.c. (Al) 



A/=i = f2e ab3 M l5 [A fl ,Xl] + T» + f 3 e ab3 m l5 [A„Xl]^ 

+^N[(X a L A,Xl - XlA.Xl) - (X a R A,X R - X R A,X a R )}T» 

+ ^{N[3X 3 A»(X l L T^ + X|T 2 ) + 3(XlA»X 3 X + X 2 L A»X 3 T 2 ) 
-2{X\A»X\ + XlA^Xl - 2X 3 L A^X 3 L )T 3 } - (L <-> R)} + h.c. (A2) 



Cti=2 = f^T'Ni^XlA^Xl + X%A lt X b R ]T^ 

+f 5 e ab3 Ni l5 {X a L A,X 3 + X 3 A,X a L + (L <- R))T£ 

+g 5 l ab N[A^X a _] + T b ^ + g 6 l ab N[A ll ,X a _]^ + h.c. , (A3) 

where the terms containing fi and g^ start off with one- and two-pion vertices, respec- 
tively. In the heavy baryon expansion, the terms containing the /j start to contribute at 
0(l/m N ). The leading order term vanishes since P + ■ 275 • P + — 0. Since we work only to 
lowest order in the 1 / m N expansion, we obtain no contribution from the terms containing 
the fi. 

For the pv 7rAA effective Lagrangians we have 

Alo=JoTM M 7^, (A4) 



£#=1 = ^fYTiTr(A^Xl) - ^Tri 5 T,Tr(A,X 3 ) 

h l - h 2 

^f n TXiTi ^f^{3T 3 (X 1 _T 1 + X 2 T 2 ) + 3(T 1 Xl + T 2 X 2 )T 3 

2a/2 2y2 

-2(f l X 3 _T l + f 2 X 3 _T 2 - 2T 3 X 3 T 3 )} + j 2 {3[(f VT 1 + f l ^ ,x T 3 )Tr{A lx X\) 
+ (f 3 7 ^T 2 + f 2 YT 3 )Tr(A tl X 2 + )] - 2(f l ^T 1 + f VT 2 - 2T 3 YT 3 )Tr(A^X 3 + )} 
+k 2 {3[(f V75T 1 + f 1 7^75T 3 )Tr(A M Xi) + (T 3 7 ^ 5 T 2 + f V 75 T 3 )Tr(/i M X 2 )] 

-2(T Vts^ 1 + f V75T 2 - 2f 3 7^7 5 r 3 )Tr(A M X!)} 
+j 3 {TV[A„ ^]+T 3 + f V L4 M , X a + ] + T a } 
+ j4 {TY[A M ,I^T 3 - TVK,i;]-T a } 
+A; 3 {f V7sK, X°] + T 3 + f Vts^, ^] + T a } 
+fc 4 {TY 75 [A /1 ,I!]_T 3 - f V 75 [A M ,X^]_T a } , (A5) 
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£a/= 2 = ja ab f a rA„T b + J& l ab T[X a R A,X b R + XlA^-fTi 

+k 5 l ab T[X a R A IM X b R - XlA.Xl^T, 

+he ab3 [T 3 i^X b + T a + T a i^X b + T 3 } , (A6) 

where we have suppressed the Lorentz indices of the A field, i.e., T v ■ ■ - T v . The vertices 
with hi, h,A contain two pions. All other vertices contain one pion when expanded to the 
leading order. At first sight the leading order term with kg in (|A6|) has no pions. However 
such a term cancels its hermitian conjugate exactly. The constants h l nA are the PV 7rAA 
Yukawa coupling constants. 

In Section 2, the leading terms in the above Lagrangians were expressed in terms of 
effective tttcNA and 7rAA coupling constants. These constants may be expressed in terms 
of the fi, gi, hi, ji and h z 7rA as follows: 

h A = ~ 2 9l + 2#2 - 03 - 3# 4 - -5-5 + -9e 



, pA+7r°7r° 



I pA+7T+7r 

h A 



9 



= 2flfi + g 3 + Q94 + 7^95 
= -^(692 + 9<74 + 2<7e 
(-65-1 - 9^4 + Ag 5 + 6# e 



a/6, 



^pAOtt+tt _ __J^ 6 ^ + _ 3^ + + 25 , 5 + 2g & 

hP AW = _^!(_ 6 ^ 1 + U g 2 + 3g 3 + 1% 4 - 2g 5 - 8^ 



^pA 1T + TT^ 



nA + + 7T 7T 



V2 
3 

3 



(65-2 - 9# 4 + 2p e 
(65-2 - 9# 4 + 2& 



— 9~( Qgi + 692 + 393 ~~ 95,4 + 25,5 + 2fl,f 



tnA T ?r 7T 



v/3, 

-g-(-%l + 12 #2 - 3^3 - 18^4 - 



806 



^/ W = -^(-60 2 + 9,4-2, 6 
= _^(_ 6 ^ + 6g 2 + 4g 5 - Ag e 



h' 



nA _ 7r+7r° 



2^1 - 25-2 + #3 + 3£f 4 + -# 5 - -5-6 
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h nA-^+ = _ 2gi _ gs _ %4 _ 2^ (A?) 



h A = K A + h 



2 

A 



^ ++A+ = -=(j + -j 6 ) - - -^-0'3 + H 



u A+A° v 7 ^ , . 4 2a/2 . 
V = -3-U0 + gJe) ^-35 

hf A ~ = ^=(jo + \h) + 2v/6j 2 + + h) + ^= (A8) 

It is interesting to note there is only one independent PV Yukawa coupling constant h& 
for 7rAA interactions. 



APPENDIX B: VANISHING LOOP CONTRIBUTIONS 

As noted in Section 3, a large number of graphs which nominally contribute to h FFF 
actually vanish up to 0(1/ A^). Here, we summarize the the reasons why. 

Consider first the corrections due to the PV vector nNN vertices. For FIG. 1 (b) we 
have 



(Bl) 

where we have used v ■ q ~ 0(l/m N ). Since we are working to leading order in the l/m N 
expansion, this amplitude does not contribute. The PV vector interactions also appear 
in Figs. I(jl,j2). The corresponding amplitude is 

tM w+W - -z-^r (h v + 2h v - -h v ) J {2n)D v . to . (k + q)(k* - ml) ~ ° (B2) 

This integral vanishes because it is proportional to [(S-v), (S-q)] + , which vanishes because 
S-v = 0. All other possible insertions of PV vector ttNN vertices vanish for similar reasons 
as either ([Bl]) or (|B2|) . In what follows, we refer only to insertions involving the PV 7rNN 



Yukawa and ttttNN axial couplings. 

The propagator corrections in FIG. 1 (gl)-(h2) vanish after integration since their 
amplitude of (gl)-g2) goes as 

W|M4=° • < B3 > 
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while the amplitude of (hl)-(h2) goes as 



The amplitude of FIG. 1 (il)-(i4) contains a vanishing integral 

~ K J (L)°v-k(k*-ml) =0 ■ (B5) 

Figs. 1 (jl)-(j2) do not contribute for the PV Yukawa coupling h w due to charge 
conservation. The remaining non-zero diagrams are Figs. 1 (a)-(f2) where the insertions 
in loops are of the Yukawa or axial interations. Figs, (fl), (f2) arises from the insertion of 
the counter terms of mass and wave function renormalization. Fig. 1 (el)-(e2) and Fig. 
2 (cl)-(c2) contribute to the wave function renormalization in Eq. fl3"B|). 

Due to the heavy baryon projection P + ■ 275 • P + = the one pion PV 7rNA vertex 
does not contribute in the leading order of heavy baryon expansion. Hence, the chiral 
loop corrections from FIG. 2 (dl)-(g4) are of higher order. Fig. 2 (hl)-(j2) vanishes 
after integration for reasons similar to (|B2]). The remaining, non- vanishing diagrams are 
discussed explicitly in Section 3. 

As pointed out in Section 2, both PC and PV two-derivative operators which conserve 
CP do not contribute to renormalization. For example, there exists one CP-conserving, 
PV such operator: 

^-Na^iD^A,, - D^N . (B6) 

After expansion, the leading term starts with three pions. It contributes via Figure 1 (a), 
at the order of 1/A X F%. Moreover the loop integration yields a factor g^ v and leads to 
zero after contraction with a^ v ' . 

Another possibility comes from insertions of PC two-derivative nucleon pion operators. 
There are three PC operators which conserve CP: 

j-Ni lb D,A»N , (B7) 



1 _-, 
A 



NA^ApN , (B8) 



^Na^A^ A^N . (B9) 

Note the first two operators are symmetric in the Lorentz indices. Only the last one arises 
from the antisymmetric operators listed in Eq. The first one starts off with one 

pion. The relevant Feynman diagrams are Figure 1 (cl)-(c2), where the PV vertex is 
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associated with h l A . Note these diagrams do not contribute at leading order of HBChPT 
due to the presence of the £75. The remaining two operators start off with two pions. The 
relevant diagrams are Figure 1 (dl)-(d2). After integration the contribution of the third 
operator reads 

~ h n e^v a Spv^ml In m n /A x F* . (BIO) 

So its contribution is zero. In contrast the second operator yields 

K{vq)ml\nmJ{K x Fl) . (Bll) 

Note v ■ q ~ l/m N . So its contribution is of order l/(A x m N ). In short, none of the 
two-derivative operators contribute to the renormalization of h n at the order to which we 
work. 



APPENDIX C: BAG MODEL INTEGRALS 

Here, we show that the four-quark bag model integrals relevant to the calculation of 
the DDH/FCDH parameters rj and y have a Taylor expansion in light quark mass around 
m q = 0. We write a bag model quark wavefunction as P8|,|T8 



^ x)= {7^)lr X ) e ^- lEt) ' (C1) 

where \ denotes a two-component Pauli spinor and where wave function normalization 
yields 

' d 3 r(u(r) 2 + £(r) 2 ) = 1 , (C2) 

where the the radial integration runs from to the bag radius, R. The four quark matrix 
elements of interest here can depend three different integrals: 



d A ru{rf , / d A rt{rf , / d 6 r w(r) 2 l(r) 2 . (C3) 
The quark radial wave functions are 

u(r) = Nj {?£) (C4) 



1/2 

uj n + m a R ) R 



f \ at I u n -m q R\ . p n r 

(r) = —N ■ Ji(-fr) , (C5) 



where 



tanp n = Pn (n=l,2,...) (C6) 

uj n + m q R - 1 
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Pn = Jul - mlR* (C7) 



N 



Pn ( Q 



\ R 3 (2luI - 2uo n + m q R) sin 2 p n 

R ~ 4ttB (C9) 

B is the bag constant and Z is a phenomenological parameter involved with the center 
of mass motion of the bag. 

For light quarks and lowest eigenmode 

ujq « (2.043 + 0A93m q R) (CIO) 
2.27/v / 4ZR3 . (Cll) 



It is straightforward to show that the bag model integrals in Eq. ( |C3[ ) have a Taylor 
expansion about m q = 0. The argument proceeds by noting that the quantities N, R, p n , 
u n and the argument of the spherical Bessel functions all have Taylor series in m q about 
m q = 0. The existence of this expansion can be seen be an explicit, iterative construction. 
First, expand u n and R: 



Jn=Y, Un,km k q (C12) 
n=0 

oo 

R=Y J R k m k q . (C13) 



n=0 



Now let m q = in Eqs. ( p6|JC7|) . Doing so eliminates all dependence on R and determines 
u n fl. Next, set m q = in Eqs. ( |C5| , |C5| ) with u n — > cj nj0 . Doing so determines i?o- Now 
expand Eqs. ( |C6| , |C^ ) to first order in m q . This step yields u n> i in terms of u) n fl and 
Rq. Expanding Eqs. (CS, |(J9[ ) to first order in m q then determines R\ in terms of k; n ,o> 



cu ni i, and i? an d so forth. Note that at any step of the recursion, the argument of any 
transcendental function is u n fl. Hence, at any order, a solution for the uj n ^ and Rk exists. 

The expansion of the bag model integrals continues by computing their derivatives 
with respect to m q and using the expansions of N, R, etc. in terms of m q as constructed 
above. Taking n derivatives of one of the integrals in Eq. ( |C3j ) yields new intregrals 
involving powers oir/R times products of the Bessel functions and their derivatives. Using 
the standard Bessel function recursion relations, the derivatives of the jk can always be 
expressed in terms of other spherical Bessel functions. Since the jk and their derivatives 
are finite at the origin, and since the radial bag integration is bounded above by R, the 
nth derivative of any of the integrals in Eq. (|C3|) is finite. Thus, each of the integrals in 
Eq. (|C3|) can be expanded in a Taylor series about m q = 0. 
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Figure Captions 



Figure 1. Meson-nucleon intermediate state contributions to the PV irNN vertex h n . 
The shaded circle denotes the PV vertex. The solid and dashed lines correspond to the 

nucleon and pion respectively. 



Figure 2. The chiral corrections from A intermediate state, which is denoted by the 

double line. 



Figure 3. Diagrammatic representation of the SU(6) W components of 
{B'M\H^(AS = 0, Figs. 3a-c correspond, respectively, to b t>v , c v , and a t:V . The 

wavy line denotes the action of . 



Figure 4. Chiral corrections to the B — > B'M nonleptonic weak decay. 



Figure 5. Quark line diagrams for the renormalization of h w due to the axial PV nirNN 
interaction. As in Fig. 3, the wavey line denotes the action of '. Fig. 5a shows a 
typical contribution to h\. Figs 5b, c denote the corresponding loop corrections to h n . 
Fig. 5(b) contains the disconnected qq insertions, while Fig. 5(c) gives a Z-graph 

contribution. 
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